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This review provides a gentle introduction to one-way quantum computing in distributed 
architectures. One-way quantum computation shows significant promise as a compu- 
tational model for distributed systems, particularly those architectures which rely on 
probabilistic entangling operations. We review the theoretical underpinnings of one-way 
quantum computation and discuss the practical issues related to exploiting the one-way 
model in distributed architectures. 

1. Introduction 

The seminal work of Raussendorf and Briegeff proposed the fundamentals of one- 
way quantum co mputin g, and these ideas were developed by follow-up papers from 
numerous sourceJ^EEl One-way quantum computing is so-called because it is driven 
by irreversible measurements performed on a large entangled state, rather than 
by reversible gates. Although the gate-based model of quantum computing uses 
entangling operations, these arc spread throughout the computation, whereas in 
one- way quantum computing all the entanglement required is present in the initial 



Many recent architectures proposed for distributed quantum computing suf- 
fer from an inability to dctcrministically entangle qubits. If entangling gates be- 
tween single qubit nodes are projective and probabilistic then a naive application 
of the gate-based approach will be prohibitively inefficient. However, projective 
and probabilistic entangling operations can be employe d quite nat urally to grow 
an entanglement resource required for the one-way 
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Conveying an 

understanding of how one grows an entanglement resource is our core goal here. 

The entanglement resources of one-way quantum computing are called graph, or 
cluster, states and our starting point is to explain what these entanglement resources 
are. We introduce two invaluable tools for describing graph states. Firstly, we out- 
line the stabilizer formalism, which uses operators to describe quantum states. Sec- 
ondly, wc introduce a graphical representation from which graph states inherited 
their name. Having described graph states, their information processing capabilities 
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will be introduced by analogy with telcportation. Graph states will be shown to be 
capable of universal quantum computing in a deterministic fashion. The determin- 
ism of one-way quantum computing will be a point of interest as it is this feature 
that was not obvious prior to the work of Raussendorf and Briegcl. Indeed, since 
the measurements on graph states necessary to perform computation have random 
outcomes it is counter-intuitive that this evolution can be harnessed to detcrminis- 
tically process information. As we shall see, certain measurements on graph states 
can be classically simulated, and so can be removed to yield a smaller, more efficient 
resource state for a given algorithm. 

Having established the basics of one-way quantum computing, we move on to 
describe how it can be physically embodied in a distributed quantum computer. 
An architecture for quantum computing is distributed whenever it has a natural 
partitioning into nodes, small quantum registers, that are interconnected by some 
network that allows communication between nodes. Each node is a quantum com- 
puter of limited size, where quantum control is excellent between that small n umber 
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of qubits. This is an apt description of many next generation proposal d 1 where 
an initial design has been found to work well on a small scale, but scaling up 
the size of a node is either difficult or simply impossible. Rather than abandoning 
an excellent small scale technology, the distributed paradigm prescribes that we 
chain up these technologies into a network of nodes. The approach is proving at- 
tractive to experimental groups because separation of nodes allows both for more 
efficient cooling and individual addressability, which become possible by virtue of 
spatial separation. In solid states systems with densely packed qubits addressing 
individual qubits remain a significant obstacle making the distributed paradigm an 
appealing alternative. 

The concept of a distributed quantum computer (DQC) is quite general, but it 
is especially suited to so called hybrid matter-optical systems. In hybrid systems 
some collection of matter qubits fo rm a node, and the inter-node e ntanglement is 
generated by an optical mechanisn ^ l 15 l 16 l 17 l 18 l 19 l 20 l 21 l 22 l 23 J2l. Such schemes 
have been experime ntally demonstrated between: ions held in separate electromag- 
netic trapJ^ESEZl anc i clouds of cold caesium atoms in separate containers^!. We 
will see that when photon detectors register a particular set of detector "clicks", 
this can herald the creation of entanglement if the photon source is unknown. Typ- 
ically, information regarding a photon's path to the detector is erased by use of a 
beam splitter, so that it cannot be determined which node was the emitter. Our 
introduction to hybrid systems is split into a two parts: firstly, a review of the 
underlying physics of atom-light interactions and linear optical transformations; 
secondly we outline a number of techniques for dealing with practical issues which 
arise in such systems. Since photons are easily lost, for example by absorption or 
detector inefficiency, robustness against this problem will play a prominent role in 
our review of these protocols. 

One-way quantum computing and distributed quantum computing are partially 
independent topics. The former is a model of computation, and the latter a physical 
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design. However, they fit together neatly to offer a complete approach to quantum 
computing. On one hand DQC gives us a method of producing entanglement based 
on photon measurements that are inherently probabilistic and projective. On the 
other hand, one-way quantum computing requires an entanglement resource that is 
simple to describe. Fortunately, several strategies exist that ensure these imperfect 
entangling operations can efficiently grow the desired graph states, and we out- 
line these in section HI We close with a brief discussion of how to tolerate other 
imperfections using fault tolerance techniques and entanglement distillation. 



2. One-way computation 
2.1. The stabilizer formalism 

The stabilizer formalisrrPSl jg a powerful tool for describing certain quantum states, 
and achieves this by a compact description in terms of operators instead of wave- 
functions. The stabilizer formalism was originally developed in the context of quan- 
tum error correction, with most of the original literature emphasizing this appli- 
cation. Our presentation of the stabilizer formalism is stripped of these details, 
providing just the essentials for understanding one-way quantum computing. 

For an operator description of a quantum state, it is sufficient to have: (?) a 
non- degenerate operator, Q, for which \ip) is a eigenfunction; and (ii) a number, qi, 
that specifies the correct eigenvalue, Q\i/>) — C[i\ip)- Non-degeneracy is the property 
of an operator that has only one eigenfunction for each eigenvalue. Hence, a single 
such operator is sufficient to identify a unique state. For simplicity, the eigenvalue 
can be absorbed into the operator to make a new operator S = (l/qi)Q. For now, 
we will call S a stabilizing operator of \ip) , similarly for any other operator satisfying 
S \ijj) = \ip). Later we further restrict the class of stabilizing operators. 

Unfortunately, as with the wavefunction, this non-degenerate stabilizing opera- 
tor may not have a simple decomposition. Instead, consider what can be deduced 
about \tp) from knowing that it is stabilized by a degenerate operator S, with j 
orthogonal eigenfunctions {|^i), IV^j ••■IV'j)} with eigenvalue +f. ft is clear that 

may be any of these eigenfunctions, or indeed any complex linear combination 
of them; that is, \ip) is within the j-dimensional subspace, J#s, stabilized by S. 
Given a set of stabilizing operators M. = {Si, S2, • ■■S'n}, IV 1 ) must be within the 
subspace stabilized by all these operators. Therefore, given enough stabilizing oper- 
ators, even if each is degenerate, we can reduce the number of possible states to one 
definite state. This set of stabilizing operators, M., forms the stabilizer of \4>)- Typ- 
ically, we consider only stabilizing operators chosen from tensor products of Pauli 
operators (X, Y and Z), the Pauli group. By restricting stabilizing operators to 
Pauli operators we limit ourselves to being able to describe a subset of all possible 
states, which we call stabilizer states. In exchange for this limitation, the stabilizer 
formalism gains a conciseness that enables it to avoid the exponential number of 
variables required to describe a general quantum state. 
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2.1.1. Commutation of stabilizing operators 

Any two stabilizing operators of a state must commute on the state-space they 
stabilize, since 

[s,-(V0, Sfcty)] \i>) = SjWjSkW) W) - s k y>)SiW IV>) (i) 

= S j (if>)\i>)-S k ^)\^) (2) 

= \1>) - W = 0. (3) 

Further, for stabilizing operations in the Pauli group this implies that the op- 
erators themselves must commute. 



2.1.2. Products of stabilizing operators 

For any two stabilizing operators, S k and Si, their product must also be in the 
stabilizer of the same state, since: 

(S k S e ) |V) = S k {St |V» (4) 
= S k |V) (5) 

= m- (6) 

As a consequence of this, in order to track the evolution of a quantum state, it 
is not necessary to individually track the evolution of all operators in the stabi- 
lizer. Instead it is possible to track only a subset of stabilizing operators, which by 
multiplication generates the rest of the stabilizer. Any minimal collection of opera- 
tors that suffices to describe the stabilizer are called generators of the stabilizer. It 
has been showrP^ that n generators are sufficient to describe all n-qubit stabilizer 
states. 



2.1.3. Unitary evolution 

Closed quantum systems, which are not interacting with an environment, are 
governed by the Schrodingcr equation. After at time A of evolving under the 
Schrodingcr equation, an initial state \tp(t)) evolves unitarily to \ip(t + A)) = 
U \ip(t)) where U is a unitary operator satisfying Wll = 11. If the initial state 
was stabilized by S, then: 

m+A)) = um)) (7) 
= us\m) (8) 

= (USU^)Um)) (9) 
= (USrf)\i/>(t + A)) (10) 

and so the new state \ip(t + A)) = U \ip(t)) is stabilized by S' — USW . In general, 
U will not map a Pauli operator to a Pauli operator, and may evolve an initial sta- 
bilizer state into some state which does not admit an efficient stabilizer description. 
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However, some unitaries will always map Pauli operators to Pauli operators. Uni- 
tary operators that satisfy this property are said to belong to the Clifford group, C, 
where the group structure entails that if C\ € C and C2 G C then the product is also 
in the Clifford group C1C2 € C. Consequently, we can describe the whole Clifford 
group by describing a subset of operators that generate the rest of the group. One 
collection of generators are the single qubit Hadamard, square-root phase flip and 
two-qubit control-phase gate: 



where in the definition of the Hadamard, H, we have used the shorthand 
|±) = (|0) ± |l))v2. The Clifford group is an important subgroup of unitary oper- 
ations, containing gates necessary to demonstrate many of the important features 
of quantum mechanics, including non-commuting operators, manipulation of quan- 
tum superpositions and creation of entanglement. One of the most striking features 
of this group of operations, however, is that unlike general unitary operations evo- 
lution under Clifford group operations often r emains classically simulable, a result 
known as the Gottesman-Knill theorerrP2E2] 



2.1A. Measurements 

What measurements can we encompass within the stabilizer formalism? If we mea- 
sure an observable, M, and get outcome rrii, then the measured system is projected 
into a state stabilized by rri^ M. Hence, we expect that measurements of Pauli 
operators can be modeled within the stabilizer formalism. Clearly, m~ 1 M is added 
to the stabilizer of the state, but does the projection remove any operators from the 
stabilizer? Recall that all operators within the stabilizer must commute with each 
other, and that all Pauli operators either commute or anti-commute. It follows that 
projection will cause all operators that anti-commute with m~ 1 M to be eliminated 
from the stabilizer. In terms of stabilizer generators, there is always a choice of gen- 
erators where no more than one generator anti-commutes with the measurement 
operator. 

In addition to simulating Pauli measurement projections we can also predict the 
outcome probabilities. The are two different cases: when the measurement operator 
commutes with the whole stabilizer group, and when it does not. If it commutes with 
the stabilizer the system must already be in an eigenstate of the observable, and 
one outcome occurs with unit probability. When this is not the case, the outcomes 
+1 and —I occur with equal probability. 




|o)(+| J - + |i)H,-, 
loxoij + iiixiij, 

|0,0)<0,0| iJ + |I,0)(l,0| JJ + |0,l)(0,f| lJ + -|l,I)(f,l| 



(11) 
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2.1.5. The Gottesman-Knill theorem 

We have loosely remarked that stabilizer formalism provides a compact description 
of stabilizer states. More formally, an n-qubit stabilizer state can be described by a 
classical bit string that scales as a polynomial in n. We can arrive at this result by a 
simple counting argument: An n-qubit stabilizer state is described by n generators, 
so we require nG bits where G is the number of bits per generator. Each generator 
consists of a ± sign (1 bit) and a tensor product of n Pauli operators (2n bitt@) , so 
G = In + 1. The grand total is n(2n + 1) bits, which is a slow scaling polynomial 
much smaller than the exponential amount of information required to describe a 
general quantum state. This argument suffices to establish the effi ciency of the 
stabilizer formalism, though further reductions have been investigate d^^ . 

This result lies at the heart of the Gottesman-Knill theorem. The theorem states 
that any quantum circuit composed entirely of gates in the Clifford group and 
Pauli basis measurements, taking as input a stabilizer state, can be simulated on 
a classical computer with only polynomial overhead. Formally the theorem can be 
stated as follows: 

Theorem 1. Any quantum computer performing only: a) initialisation of qubits in 
stabilizer states, b) Clifford group gates, c) measurements of Pauli group operators, 
and d) Clifford group operations conditioned on classical bits, which may be the 
results of earlier measurements, can be perfectly simulated in polynomial time on 
a probabilistic classical computer. 

On the basis of the preceding sections, this should seem plausible to the reader. 
One way to see that the Gottesman-Knill theorem must hold is to consider the 
effect of such a circuit on the stabilizer of the initial quantum state. This initial 
state has by definition an efficient stabilizer description. In general, measurements 
can be postponed until after the unitary part of a quantum circuit, by adding a 
single ancillary qubit per measurement. 

From section T2. 1.31 the stabilizing operators, {5fc}, of the initial quantum state 
evolve under an operator U as USk . If U is an element of the Clifford group 
and {Sk} are Pauli operators, then the resultant generators will also be within 
the Pauli group. Thus the quantum state produced as a result of the circuit must 
also be a stabilizer state. As the size and form of the set of generators of the 
stabilizer describing the quantum state remains constant, the number of parameters 
describing the state also remain constant. This is in stark contrast to the general 
case, where the number of free parameters describing a quantum state can grow 
exponentially with the number of quantum gates. For a formal proof of the theorem 
and exact compu tational reso urces required for such a simulation we refer to the 
relevant 

literature^™™ 



a The factor of 2 arises because we use 2 bits to describe each Pauli operator, e.g. ]1 = (0,0), 

x = (o,i), y = (i,i),z = (i,o). 
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2.2. Graph states 

2.2.1. The constructive definition 

In mathematics, a graph is a particular kind of diagram composed of points (graph 
vertices) connected by lines (graph edges), where some simple examples are shown 
in tables [1] and [2l In quantum information graph states are a particular class of 
quantum states that can be represented by a graph. There are numerous definitions 
of a graph state, but the simplest definition to grasp is the constructive definition, 
where a graph specifies a procedure for constructing the corresponding quantum 
state. Firstly, each graph vertex, i, represents a qubit prepared in the state !+},•• 
Secondly, each graph edge between two vertices, i and j, represents the application of 
a two-qubit controlled-phase gate, CZf (defined in Eqn. (fTTjl ). All CZ operations 
commute, and CZ is unchanged on interchange of the control and target qubit. 
Therefore, the edges do not require a time ordering or bias toward one qubit. 

Cluster states are a specific class of graph state for which the corresponding 
graph is a regular square lattice. These were initially introduced by Briegel and 
Raussendorf in the context of studying many-body entanglement!^, but as we shall 
see in section f2. 31 also provide a universal resource for measurement based compu- 
tation. Although the constructive definition specifies a unique quantum state for 
each graph, the corresponding wavefunction rapidly becomes cumbersome as the 
graph grows in size. 

2.2.2. Graph states as stabilizer states 

As we have seen in section 12.1.51 any stabilizer state to which a Clifford group 
operator is applied remains a stabilizer state. The constructive definition defines a 
graph state as the result of just such a process, and so graph states are stabilizer 
states. As the initial states to which the CZ operators are applied are stabi- 

lized by Xi, by evolving the generators of the stabilizer as discussed in section [2.1.31 
we obtain the corresponding operators for the graph state. This yields an alternate 
definition for the graph state corresponding to some graph G as the state for which 
the set of operators 

s=x l n z i (12) 

jeN G (i) 

generate the stabilizer. Here Nc(i) denotes the set of vertices connected by an edge 
to vertex i, which is called the neighborhood of vertex i. 

So, all graph states are stabilizer states. What about the converse statement: 
are all stabilizer states also graph states? Consider the two-qubit state stabilized 
by X1X2 and Z1Z2] this is an example of a stabilizer state that cannot be gen- 
erated by the constructive definition, and so the converse statement is not true. 
However, it can be transformed into a constructively defined two-qubit linear chain 
by the application, to either qubit, of a Hadamard gate. Since the Hadamard is 
a local Clifford operation, there exists a graph state with the same entanglement 
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Table 1. Three representations of graph states for: (a) two-qubit linear chain; (b) 4-qubit linear 
chain; (c) a 5-qubit non-trivial graph. 



Graph 

representation 


Wavefunction 
representation 


Stabilizer 
generators 


(a) 

1 2 


\%) = czl\++) 

= (l/V2)(|0+) + |l-» 


X%Z\ 


• • 


(b) 

12 3 4 


\%) = cz\czlcz\ |+ + ++) 

= (1/V2)( |+00+) + I+01-) 
I-10+) + 1-11-)) 


X2ZIZ3, 
X3Z2Z4, 
X4Z3 


• m m • 


(c) 

12 3 4 


|sf c ) = cz 2 i czlczlcz'!\+ + + + +) 

= (1/V2) (|+00 + +) + |+01 -+) 
+ |-10 + -) + |-ll--)) 


x x z 2 , 

X2Z1Z3Z5, 
X3Z2Z4, 
X4Z3, 
X5Z2 


* 1 * * 

l 5 



properties as our example stabilizer state. Surprisingly, this is always the case: all 
stabilizer states are local Clifford equivalent to some graph state^EH. A loose way 
of understanding this result is to consider how the constructive definition provides 
the CZ\ gates, and so only the local gates Hi and Pi are required to complete the 
generators Clifford group. However, this mapping is not necessarily one-to-one, as 
shown in table [2j a single stabilizer state may be equivalent up to local unitary 
operations to many different graphs. 

In table [2] the required local unitaries have been explicitly added so that the 
different representations arc strictly identical. However, if we are interested in the 
required entanglement resource then local unitaries arc unimportant, and are com- 
paratively easy to implement. Therefore, different topology graphs can be sorted 
into local unitary (LU) equivalence classes. Hein et al have characterized all of the 
equivalence classes for graphs with up to 7 vertices^. They found that out of 995 
graphs there are only 45 classes of LU-equivalent graphs. Local Clifford (LC) equiv- 
alence is a special case of LU-equivalence, in which equi val ence is only considered up 
to operations in the Clifford group. Van den Nest et showed that the problem 
of identifying whether two graphs arc LC-cquivalcnt is reducible to a purely graphi- 
cal problem; two graphs are LC-equivalcnt if and only if there exists a series of local 
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Table 2. Example stabilizer states with different equivalent graph representa- 
tions. Vertices labeled with operators are equivalent to the stabilizer state when 
these local unitaries are applied to the vertex after having followed the constructive 
prescription. 



Equivalent graph 
representations 


Stabilizer 
generators 


f\ 12 - 12 
(a) H . . - . . H 


Z1Z2 


2 j>HPH 

w — 3 - P Ap» 

1 3 


X1Z2, 
X2Z1Z3, 
X3Z2 



complementations that transforms one graph into the other. A local complementa- 
tion is an operation that can be applied to any vertex i. When applied to vertex i, 
every neighbour of i gains a graph edge with every other neighbour. However, since 
(CZ*) 2 = 1, if the neighbours were already connected, the two graph edges cancel 
each other. In table [2Jd we transform between the two graphs by performing a local 
complementation on vertex 2. As vertex 2 has 1 and 3 as neighbouring vertices, 
the action of complementation connects or disconnects 1 and 3. This leaves the 
problem of finding an algorithm to work through all possible combinations of local 
complementations, so as to determine whether two graphs are L C-e quivalent. Van 
den Nest et showed that an algorithm developed by BouchelP^ — initially for 
purely theoretical purposes — will efficiently solve this problem in polynomial time. 
This algorithm proves extremely useful when planning how best to implement an 
algorithm using graph states. Many architectures do not allow for many gates to 
be simultaneously executed on a given qubit, meaning that edges must be created 
sequentially, as is the case for many schemes base d on optical measurements, such 
as double-heralding^^ and repeat- until-succesJ^. In these cases it is sensible to cy- 
cle through LC-equivalent graphs to minimize the number of entangling operations 
necessary to construct the state. 

2.2.3. Pauli measurements on graph states 

When graph states undergo Pauli measurements the transformation obeys a set 
of graphical rule s, which were discovered independently by Hein et a$* and by 
SchlingemanrP^. Here we will only describe the rules, up-to local Clifford operation, 
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Table 3. Transformation rules for graph states undergoing local measurements of Pauli group 
operators. These describe the resultant state up to local Clifford operations. Here we adopt 
the convention of denoting sets by uppercase letters and elements of a set using lowercase 
letters. The equations in the table use E(A, B) to denote the set of edges between sets A 
and B, E{A, B) = {{a, b} £ E : a £ A,b £ B,a ^ b}, and the operation EAF is denned as 
EAF = (E U F) — (E n F). In the examples, vertex a is highlighted in white. 



Rule 


Example 


For a Z-basis measurement of vertex a: 
(1) Remove vertex a and edges from a. 




Formally: G' = G - {a} 




For a y-basis measurement of vertex a: 

(1) Local complement vertex a; 

(2) Remove vertex a and edges from a. 

Formally: G' = GAE(N a ,N a ) - {a} 


Y 


For an X-basis measurement of vertex a, first choose a 
vertex b (grey) that is a neighbor of a, and then: 

(1) Complement the neighbors of b with the neighbors of a; 

(2) Complement the mutual neighbors of vertices a and 6; 

(3) Complement vertex b with the neighbors (except b) of a; 

(4) Remove vertex a and edges from a. 

Note that, depending on which vertex is chosen as b, the 
final graph will appear different. However, we are only 
concerned with graphs up-to local unitaries, and different 
choices of the vertex b lead to final graphs differing by 
only local unitaries. 


m 


Formally: O = GA E(N b , N a )AE(N b HN a ,N b n N a ) 
AE({b},N a -{b})-{a} 
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Fig. 1. a) Cluster state fragment implementing a Toffoli gate. Qubits measured in the computa- 
tional basis measurements have been removed, b) The minimal graph state obtained by removing 
the Clifford part of the first graph. The unmeasured qubits retain their index from (a). Remov- 
ing the Clifford part of an algorithm in this way can dramatically reduce the number of qubits 
necessary to implement the computation. 



which will mean that the outcome of the measurement is unimportant. The rule 
for Z measurements is rather simple, the measurement removes the corresponding 
vertices from the graph. Since X and Y Pauli measurements differ from Z by a 
local Clifford operation, we expect they begin with a set of local complementations 
before removing a vertex from the graph. Given a graph 67, the rules by which the 
structure of the graph is changed by local Pauli basis measurements of vertex a is 
summarized in table 12. 2T3l where we use of the notation of Hein et aPUo formalize 
the transformation. 

Our motivation for studying graph states is that measurements on them en- 
able quantum computation, which we describe in the next section. Since the effect 
of Pauli measurements can be easily simulated, it is interesting to ask what this 
might mean in terms of the underlying computation being performed. Given the 
duality between Pauli measurements and Clifford group gates, the reader may not 
be surprised that Pauli measurements on a graph state execute the Clifford gates 
within a computation. The net result of this is that it is possible to remove the 
entire Clifford part of an algorithm, transitioning from a regular cluster state, to a 
more algorithm specific, less regular graph state. A graph fro m w hich the Clifford 
part has been removed is referred to as a minimal graph statS^, and an example 
of this removal process is shown in Fig. [TJ This ability to remove any Clifford part 
of an algorithm is strong motivation for considering general graphs, rather than 
remaining restricted to cluster states. 

2.3. One-way quantum computation 

In 2001, Rausscndorf and Briegcl proposed an innovative way to exploit the 
evolution o f clu ster states under measurement to perform universal quantum 
computation^^. They proposed measurement patterns to perform arbitrary sin- 
gle qubit unitaries, and CNOTs between logical qubits. These patterns can then be 
superimposed on a regular square lattice cluster state. Qubits not used are simply 
measured out in the Z basis, effectively removing them from the graph. 
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2.3.1. Quantum teleportation 

Quantum teleportationSSl is not only an interesting application of quantum com- 
munication, but as we will see here, also offers a useful picture for understanding 
how quantum computation can be accomplished through carefully chosen measure- 
ments on a sufficiently large cluster state. While not as powerful as the stabilizer 
formalism described earlier, this approach is particularly intuitive and is useful for 
gaining some feeling for how these computations progress. A detailed explanation 
of this approach can be found in Refs. El [4TJ 



1+) 
1+) 



IT 



mi 

m 2 

Z mi X m2 |^ 



Fig. 2. Quantum teleportation of a single qubit. A quantum circuit used to prepare the entangled 
state is boxed in grey. The top two qubits are held by one party, say Alice, and the bottom qubit 
is held by a second party, say Bob. Initially, Alice holds \ip) and without using any quantum gates 
that interact with Bob, is able to transfer it to Bob up to some Pauli corrections. Alice needs to 
transmit two classical bits of information for Bob to be able to apply the Pauli corrections and 
retrieve the teleported state. 



Let us start by examining the teleportation protocol, which is best depicted as 
a circuit diagram in Fig. [21 The steps are as follows: 

(1) A maximally entangled pair of qubits is generated, and shared between two 
parties. 

(2) The qubit to be transfered interacts with one half of the entangled pair by 
means of a controlled phase gate. 

(3) Both qubits held by the first party are then measured in the X basis. 

(4) The measurement results are communicated to the second party classically. 

(5) The second party performs a Pauli operator on their qubit dependent on the 
previous measurement results in order to obtain the teleported state. 

The quantum teleportation protocol can be seen as a concatenation of two 
iterations of the two-qubit teleportation illustrated belo w in Fig. [3] This two-qubit 
teleportation procedure was proposed by Zhou et aJ32l. 



l+>- 



H 



mi 



X mi H|V>) 



Fig. 3. The two-qubit quantum teleportation protocol. 
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We can follow this two qubit protocol mathematically. First, the entangling 
gate is performed between the qubit carrying the state to be teleported, \tp) = 
a\0) + j3 11), and the second qubit prepared in the |+) state: 

CZI\tP) a \+) b = a\Q)\+) + p\l)\~) ■ (13) 

We can emulate the effect of the X basis measurement, by instead performing a 
Hadamard rotation which will later be followed by a Z measurement: 

B A CZ B A \i') A \+) B = ^=K|0, +) + |1, +)) + m -) - |i, -))] (14) 
= -L[|0> (a |+) + /? |-)) + |1) (a |+) - /? |-))] (15) 

= -l(|0> + |l>X)( a |+)+/?|-)). (16) 

After a Z measurement with outcome mi, the state is projected by {m\\: 

( mi \ A H A CZ B W A \+) B = X m - (a\+)+/3\-)) (17) 

= X m -H\^) B . (18) 

Here the teleportation protocol introduces a Hadamard gate onto the transmitted 
qubit. The evolution introduced in this way is the key to performing computation 
via measurement. 



2.3.2. Single qubit gates 

In the one-way model, single qubit gates are implemented by way of a variation of 
the two-qubit teleportation protocol: As any Z rotations added to the state, \ip), 
commute with the controlled phase gate and so can be absorbed into the measure- 
ment basis. Rotating the X-measuremcnt about the Z-axis enables a qubit mea- 
surement along any axis in the X-Y plane. Fig. U illustrates this process, showing 
three equivalent quantum circuits leading to the same output state. 

By generalizing the measurement operator to allow for any X— Y plane measure- 
ment basis, it is possible to apply an arbitrary operation of the form Jg = He , 
for any real parameter 9. This is equivalent to applying an arbitrary Z rotation 
followed by a Hadamard gate. By concatenating three such teleportations, the re- 
sulting operator will by of the form: 

J^JcJe = He^ z He^ z He l9Z (19) 
= He* z e< x e wz 

Any single qubit unitary operation can be implemented in this way through suitable 
choice of 9, £ and £, a result which follows from the Euler angle decomposition of 
operators in SU (2)^. The initial entangling operations used in this protocol can be 
seen as preparing a graph state, upon which single qubit measurements are made. 
As a result, measurements on a linear segment of a graph state can be used to apply 
any single qubit unitary operation, as depicted in Fig. HJd). 
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Fig. 4. Rotation of a single qubit via measurements, a) The two qubit teleportation protocol, b) 
The rotation can be separated from the logical qubit. c) The z rotation commutes with the CZ 
gate, d) The measurement pattern to implement any single qubit gate exploiting the Eulcr angle 
formula, along with the initial state of each qubit. Here we use a diamond vertex to denote the 
input qubit and a square vertex to denote the output qubit. 



2.3.3. Entangling gates 




Fig. 5. Entangling gate between two logical qubits via measurements. 



Just as single qubit rotations about the z axis commute with the controlled phase 
gates used to construct the graph state, controlled phase gates also commute with 
each other, and so can simply be added to the graph to perform entangling gates 
between logical qubits. This yields another variant on the two-qubit teleportation 
protocol, which is shown in Fig. [5l 

This variation on the teleportation protocol is depicted in Fig. [6^. A further 
variant on this approach is to add an additional node in place of the edge used to 
implement the controlled phase gate, as shown in Fig. EJb) and (c). This allows 
freedom of choice in the location of entangling gates, as measuring the intermediate 
qubit in the Y-basis recovers the entangling gate, while measuring it in the Z- 
basis leads to an identity operation. Conveniently, as the graph segments used for 
both single qubit gates and entangling gates are planar, and form a universal set 
of operations for quantum computation, it follows that any quantum circuit can 
be implemented by simply superimposing the measurement patterns corresponding 
to individual gates on a regular square cluster state. Additional qubits present in 
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the original cluster state but not used in these gates can be removed simply by 
measurement in the if-basis, as discussed in section [2. 2. 31 



a. b. c. 




Fig. 6. a) Measurement pattern for implementing the entangling gate shown in Fig. [5] This can 
be replaced with an arrangement of three qubits which allows the entangling gate to be performed 
using the measurement pattern in (b), or can be used to implement single qubit gates only, by 
using the measurement pattern in (c). As before, we use a diamond vertices to denote the input 
qubits and square vertices to denote output qubits. 



3. Physical implementations 

We now turn our attention to how physical systems form a distributed architecture 
capable of producing entanglement between nodes. We focus on hybrid matter- 
optical systems where matter systems embody the logical qubits and photons are 
used to generate entanglement between nodes. The role of logical qubits is more 
suited to matter systems as matter storage is less problematic than photon storage. 
However, photons can travel long distances without suffering significant decoherence 
from interactions with the environment. Photons can be very robuslEI to decoher- 
ence because in free space photons do not interact and when travelling through 
materials the degree of non-linearity is typically extremely weak. However, generat- 
ing entanglement by unitary evolution requires a non-linear Hamiltonian. Altho ugh 
there are some schemes for quantum computing that utilise weak non-linearities^, 
this is not essential as non-linearities can also be induced by measurement. Al- 
though there are a broad spectrum of optical methods for producing entanglement, 
we have in mind implementations that use only linear optics and non-number re- 
solving photon detectors, as these present the lowest technological barriers. 

We begin our introduction to the physics of these systems by introducing the 
physics of free photons and linear optical elements such as beam splitters, phase 
shifters, and wave-plates. Next, we introduce atom-light interactions and how lasers 
can be used to manipulate matter qubits. 



b The degree of robustness depends on which photon degree of freedom is used to encode qubits. 
For instance time-bin, spatial mode and frequency encoding are more robust than polarization 
encoding, as waveguide materials can cause polarization to drift. 
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3.1. Photons and linear optics 



By quantization of the free space electromagnetic ( EM) fie ld within a cubic volume 

'k\ ( + I s 



of space V, it is well known that the Hamiltonian i 



.A' 



k.c 



where at - and a fc | are creation and annihilation operators for a photon of wavevec- 

tor, k = kk, and horizontal (e = h) or vertical (e = w) polarization^. The sum over 
wavevectors covers all frequencies within the quantization volume and three inde- 
pendent Cartesian directions. Note also that we take natural units where H = 1. It 
is useful to eliminate the need for the free space EM Hamiltonian by absorbing the 
time evolution into the photonic operators, so that r — > a? ?exp(— ikt/c). Each 
photon corresponds to a plane wave oscillation of the electric field accompanied 
by the orthogonal magnetic oscillation. Conversely, the electric field operator for 
a point in space r at time t is related to these photonic operators by a Fourier 
transform: 

E(r, t) = i ^- — ^ e ^a^ ?exp(— ikt/c + ik.r) — at ?exp(ikt/c — ik.r)j , 

(21) 

where £o is the permittivity of free space. The state of electromagnetic (EM) field 
containing a single plane wave excitation is simply at ? |vac) where |vac) is the 
ground (vacuum) state of the EM field. However, we will typically be dealing with 
localized single photons that arc a superposition of plane waves: 

= Y^c(k,i) a y (22) 

fe,e 

where normalization of the variable c(fc, e) ensures that |vac) is an eigenstate of 
the number operatoj^ with eigenvalue 1. A localized photon is not an eigenstate of 
the free photon Hamiltonian and hence evolves in time non-trivially. As we would 
expect this evolution results in propagation. From the electric field operator we can 
see that this observable is the same for a later time t' = t + At at the shifted spatial 
point r' = f + tk/c, which is consistent with a photon travelling across space at the 
speed of light. 

Having given a brief introduction to optics in free space, we now turn to the 
effects of linear optical elements on localized single photons. Linear optical effects 
are those where Heisenberg evolution of the creation operators is described by a 
mapping of vectors of operators a- m = (a\, 02, ...ajy)^ that is unitary: 

«L = Ud\ n . (23) 

c Note that, photon polarizations arc always orthogonal to a photon's wavevector e.k = 0. 
d The number operator counts the number of excitations and is n = J^r -.aj =■ 
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Fig. 7. Some commonly used linear optical elements and the corresponding operator transforma- 
tions. Phase-shifters and waveplates act on a single spatial mode, though the latter has a polar- 
ization dependence. Beam splitters and polarizing beam splitters map two input spatial modes, 
ai and 0,2, to two output spatial modes, 03 and 04. For brevity the polarizing beam splitter is not 
completely general, but rather a perfect reflector of vertical light. For the waveplate, the angle 
(p describes the relative dispersion of horizontal and vertical polarized light, whereas 8 describes 
the dispersion of polarized light in a different basis. For the beam splitter, <f> describes the phase 
acquired passing through the beam splitter, and cos 2 (0) is the percentage of transmitted (rather 
than reflected) light. 



Note that this criteria for linear optical unitaries excludes some more general unitary 
evolutions of the EM field. Consider, spontaneous down-conversion where a single 
photon splits into two photons of lower energy at — > a:- a- . Although spontaneous 

k h\ K2 

down-conversion is part of a unitary process, it is nonlinear as the mapping of 
photon creation operators can take one creation operator to two creation operators. 
In Fig. [71 we list some common linear optical elements and their effect on photonic 
operators. All of these linear optical transformations can be produced by using 
materials of different reflectivity and refractive index, with polarization sensitive 
transformations requiring materials with polarization sensitive properties. 

In many circumstances linear optics produce results that are intuitive from a 
classical perspective. However, one notably quantum effect is photon bunching, 
which occurs when otherwise identical photons impinge on different input ports of 
a 50/50 beam splitter (9 = 7r/4 and <\> = 0), so that: 

a\a\ I vac) — ► 7^( a 3 + a \){ a \ ~ a \) l vac ) 1 (24) 
-> 2(44 -44) I vac). 

Hence two identical photons tend to bunch together and leave a beam splitter 
from the same output port rather than following a classical distribution. In any 
actual experiment photons are never perfectly identical, so the effect is typically 
measured by the Hong-Ou-Mandel dipP^that quantifies the deviation from classical 
behaviour. 

Later we shall see that bunching effects can be used to produce entanglement 
by enabling a partial measurement in the Bell basis. For a Bell state to be well 
defined there must be two different photonic modes (e.g. different polarizations or 
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frequencies) at each beam splitter input. It is easy to verify that the singlet )) 
is the only Bell state that can lead to unbundled photons: 

■^7=( a \,h a 2,v - a i,v a l,h) l vac > -» -^( a l,h a l,v - a i,v a \,h) l vac ) • ( 25 ) 

This enables a partial Bell measurement because detecting antibunched photons 
will distinguish the singlet from the other three Bell states; unlike a full Bell state 
measurement that distinguishes all four Bell states from each other. 

Since a partial Bell measurement introduces a probabilistic element to entangle- 
ment generation, a natural question is whether there is any way of performing a full 
Bell measurement with linear optics. However, Vaidman and YorarM-^ showed that 
linear optics cannot be used to perform me asur ements in a maximally entangled ba- 
sis with probability 100%. Lutkenhaus et extended this proof to cover schemes 
that use additional ancillas, forward feeding of measurement results or number- 
resolving detectors. On the other hand, it is relatively e asy to design app aratus to 
project onto the {|00), |11), |01) + |10), |01) - |10)} basi J 49 ' 50 ' 51 ' 17 ! 18 ' 52 !. 



3.2. Matter-light interactions 

Since linear optical technology is most mature for optical wavelength photons, we 
will typically be interested in electronic transitions at optical wavelengths. Within 
a matter system there may be several different available transitions, but here we 
consider the simplest possible system where an excited electronic state |e) decays 
to a ground state |1) via the emission of a photon. If the purely EM contribution 
is removed from the Hamiltonian by absorbing the evolution into the photonic 
operators, then the remaining Hamiltonian for the system is: 

H = -ffmatter + Hint , (26) 

where -f/mattcr is simply the difference in the energy, w, between |1) and |e): 

-^matter = LU(T + <7~ , (27) 

which has been expressed in terms of the atomic raising and lowering operators, 
a + = |e)(l| and er~ = |l)(e|. The term H mi is the interaction Hamiltonian between 
the two systems. Typically the strongest matter-light interaction is the electric 
dipole transition^, where the energy arises from the interaction of the electric dipole 
moment with the electric field: 

Hi„t = - J M(0 ' W, t)dr, (28) 

where fi(r) is the electric dipole operator. Since optical wavelengths are much longer 
than atomic distances, we can make the electric dipole approximation that assumes 



c Transitions may be electric dipole disallowed, but may still occur by the weaker magnetic dipole or 
higher order electric contributions such as electric quadrupole. Nevertheless, the essential physics 
is captured by considering only electric dipole transitions. 
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the electric field is constant over the whole volume where the dipole moment is non- 
vanishing. When the matter system is centered at the origin, this gives a simplified 
expression for the electric field that factors outside the integral: 

i 

H^t =~i^2 (j~vi) (/ ^ d ^j ■ g ( a fc,f ex PHI fe l*/ c ) - al s exp(i\k\t/cf) . 

k,e 

(29) 

The electric dipole operator is simply the distance of charge from the origin: 

fi(f) = -ef, (30) 
= -e{x€\ + y&2 + zez), 

where e is the charge of an electron and ej are unit vectors in Cartesian coordinates. 

Since we are only interested in the action of J fx on the {|1) , |e)} subspace, we 
rewrite the dipole operator in terms of transition elements: 



^l(r)^b{r)r.eidf) , (31) 



K a,be{l,e} 



and hence we require some knowledge of the electronic wavefunctions ipi(f) an d 
il) e (r). In the solid state the problem can be quite involved as an atom in a crystal 
lattice will not have a rotationally symmetric Hamiltonian, and hence will not have 
rotationally symmetric energy eigenstates. 

However, the essential physics is captured by considering the simple energy 
eigenstates of an atom in free space. The charge of the atomic nucleus provides a 
rotationally symmetric Hamiltonian for the orbiting electron, whose energy eigen- 
states are charge balanced I^WI = IV'iC - This symmetry entails that the fol- 
lowing integrals vanish: 

' ^* 1 {r)i> 1 (7)rI i dr = 0-,Vi; (32) 
tl>l(r)ij} e (r)r.eidf = 0; Vz; 

so that the non-transitioning elements (|1)(1| and |e)(e|) also vanish. Further 
progress is easy if the ground state is an 5-shell electronic state and the excited state 
is a P-shell electronic state with angular momentum +1 in the £3 (z) direction^. 
For such a transition the symmetry of these states entails that the £3 transition 
elements vanish, and the £\ and €2 elements only differ by a factor of i. Therefore, 
for this simple system the electric dipole operator becomes: 

n(r)dr= (aa + + a*a~)ei + i(aa + — a*a~)e2, (33) 



Rotational symmetry can be reintroduced by also accounting for the presence of the P-shell 
electronic state with —1 angular momentum. 
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where the variable a depends on further details of the electronic wavefunctions. 
The phase of a is trivial as it follows from the global phase of the electronic wave- 
functions, tpi and ipj, and hence we take a = \a\. With this expression for the 
electric dipolc operator the Hamiltonian is still too unwieldy to admit any analytic 
solutions. However, the Hamiltonian is greatly simplified by neglecting terms that 
rotate quickly in the interaction picturqf|. We have already brought the creation 
operators into the interaction picture, and absorbing the effect of 7J ma ttor into the 
matter-system operators gives <r + — ► a + cxp iwt and o~ — > a~ cxp -1 "*. The quickly 
rotating terms, a + a l ex p l ( UJ + k / c ) t anc ] cr~a cxp~^ w+fc / c '*, have negligible impact on 
long term evolution. Neglecting them is the so-called rotating wave approximation 
and gives the interaction Hamiltonian: 

H in t = -iW k a+b k cxpHA fc t) - a~b\ cxp(iA fc t)^ , (34) 

where = (ui — k/c) is the detuning, Wk is the coupling strength: 

' |fc| 

£ Vct 



W k = 2a (-^-) , (35) 



and the new photonic operator is: 

= r { a ke 3 ,ei + i a ke 3 ,e 2 + a ke 2 ,e{ + ^feel.ei) ■ (36) 

This Hamiltonian shows that atomic excitation absorbs a photon, and that the 
reverse process of atomic relaxation always emits a photon. Note also that these 
processes conserve angular momentum, as the photonic operators arc +1 angular 
momentum eigenstatct0, so when emission occurs the photons carry away the an- 
gular momentum that is lost when the electron transitions from a P to S electronic 
state. 

In deriving this Hamiltonian we have made three assumptions: (z) The electric 
field is uniform over the matter system; (ii) rapidly oscillating terms can be ne- 
glected; (Hi) the electronic wavefunctions are free space atomic wavefunctions. The 
last assumption essentially gave us the polarization features of our Hamiltonian, 
and for matter qubits embedded in a crystal lattice we can expect a similar result 
except with different polarization features. For example, in both quantum dots and 
nitrogen- vacancy centres in diamond (NV centres) the X — Y plane rotational sym- 
metry can be broken and there are often t wo separate en ergy eigenstates that decay 
via horizontal and vertical polarized lighl ^ l ^'^^ l ffi In NV centres the crystal 



g Wherc evolution due to Hem and //matter have been absorbed in the the photonic and atomic 
operators, leaving iZi nt as the only remaining term in the Hamiltonian. 

h Angular momentum conservation is a consequence of the invariance of physics under rotations 
of the reference frame. Similarly, angular momentum eigenstates are rotationally symmetric up 
to a phase. For an eigenstate of eigenvalue m, the eigenstate acquires a phase of exp(im07r) for a 
rotation of 9 degrees about the angular momentum axis. Notice that the operators satisfy this 
rotational symmetry 
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lattice^ breaks the symmetry, whereas in quantum dots the dominant effect is the 
eiliptical shape of typical dots. 

Having derived Hamiltonian Eqn. (|34p . we now need to consider how to solve the 
Schrodinger equation for systems evolving under this Hamiltonian. In general, this 
is a difficult problem. However, by assuming that the coupling strength is approx- 
imately constant over the regi on whe re the detuning is small, Wk ~ W u , one can 
derive Weisskopf-Wigner This analytic solution gives an exponential 

decay of the excited state at a rate of , and will also predict many spectral 

features of the resulting photon. 

However many proposals for a mature hybrid technology employ the cavity 
effects of quantum electrodynamics, which occur between reflecting surfaces. In 
this context the constant coupling strength assumption is far from true. These 
cavities dramatically affect the coupling strengths of different frequency photons, 
depending on whether the frequencies are resonant with the cavity geometry. The 
cons truction o f such a cavity can take m any forms, including spherical mirror cav- 
ities | 59 | 60 | 61 1 microp illari^, microdisks^!, microtoroidal resonators^] and pho- 
tonic 

crystalfPlMlSI 

Interesting phenomena occur when a matter qubit is placed 
within a cavity, and under certain very general conditions is described by the Jaynes- 
Cummings modei^l In this model we assume that one resonant cavity frequency 
dominates evolution and that all other frequencies can be neglected. Assuming that 
this frequency is on resonance with the electronic transition, we have: 

#jc = igui{^~^ w - <y + b w ), (37) 

where the photon only extends across the volume of the cavity, and hence the 
relevant quantization volume is that of the cavity; we use g u rather than W u to 
identify variables with dependence on cavity volume. Furthermore, the photon is 
now a standing wave rather than a plane wave as the cavity walls enforce boundary 
conditions on the EM field. Consequently, the coupling strength is modified by a 
factor of sin(wa;/c), where x is the distance of the matter system from either cavity 
wait 

Important effects of this Hamiltonian include Rabi oscillations and the Purcell 
effect. Rabi oscillations occur when the system is initialized in a state with a single 
excitation of either matter qubit or cavity mode, and they correspond to oscillations 
of the excitation between the electron and cavity. When the system is not lossy, 
the Rabi oscillations complete full cycles with a frequency g w . When the system 
is lossy, the Rabi frequency increases and stimulat ed em ission from the cavity is 
enhanced (a phenomena known as the Purcell effedP™). 



'The NV defect has trigonal discrete symmetry, which breaks the continuous rotational symmetry 
of a free atom in space. However, the effects of symmetry breaking are amplified by strain inside 
the crystal lattice. 

J This assumes a ID model of the cavity, though a full 3D model will include a similar sinusoidally 
varying effect on coupling strength. 
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This introduction of the Jaynes-Cummings model has been idealized as we have 
assumed that photons do not escape from the cavity or get absorbed by the cavity 
walls. The quality factor of a cavity, Q, is a measure of that cavities' performance 
at retaining a quanta of energy. The value of Q is inversely proportional to the 
probability, per Rabi cycle, that the excitation is lost from the system. 

For the purposes of quantum computation, photons lost by absorption into the 
cavity walls are undesirable. However, leakage of photons into modes outside the 
cavity can be desirable, provided that the leakage is into one particular, chosen, 
mode. This can be achieved by thinning a region of one reflecting surface, which 
increases the probability of tunnelling via this route. 

For DQC applications we will need our cavity to leak photons into a mode that 
is monitored with photon detectors. If we monitor a joint matter-cavity system for 
the emission of photons, and none are detected, then its evolution can be described 
by a non-Hermitian conditional Hamiltoniarl-^4 

ffleakyJC = HjC — 2^^' (^®) 

where J is known as the quantum jump operator, 

J = Jnb u , (39) 

with k quantifying the leakage rate of the cavity. It is the non-Hermitivity of the 
additional term that is responsible for the irreversible evolution of the system. 
This irreversible decay also reduces the norm of the wavefunction N = (^l^) at 
a rate N that represents the probability density of detecting a photon. In general, 
any measurement event associated with a jump operator J will occur with proba- 
bility (\&| jt J\&). If a detection event occurs, the system is projected by the jump 
operator. 



3.3. Local rotations 

In this review we are interested in how entanglement can be produced between 
different nodes of a distributed quantum computer. However, this presupposes that 
we already have excellent control over the quantum systems within a node. Indeed, 
we will regularly assume that we can perform arbitrary local operations to high 
fidelity, with imperfections in local operations analysed as small noise contributions. 

Different physical systems allow different mechanisms for performing local uni- 
taries. For example, when we wish to rotate between a ground state |1) and an 
optically accessible excited state |e), this can be accomplished by applying a laser 
at a frequency resonant with the transition frequency. The angle 9 of the rotation, 

k Here we use irreversible to mean not cyclic, in contrast to Rabi oscillations. When the Jaynes- 
Cummings approximation is not made, we observe this kind of irreversibility in free-space photon 
emission. In the free-space case, irreversibility occurs because we couple to an infinite number of 
field modes, and hence any oscillation would take an infinite amount of time. Since, the Jaynes- 
Cummings model involves a finite number of modes, this irreversibility must be added in by hand. 
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e.g. |1) — > cos(8) |1) + sin(#) |e), is determined by the laser amplitude integrated 
across the duration of the pulse. Of particular interest are so-called ir pulses that 
pump |1) — > \e). Note also, that if each node contains more than a single qubit, 
then arbitrary rotations between those qubits are also considered to be local. For 
example, a control NOT gate between an electron spin a nd a ne ar-by nuclear spin 
is local to the node and can also be implemented opticalljE^ESl. In contrast, gates 
between qubits in different nodes can not be implemented so easily, and require an 
entanglement generation protocol. 

Since we assume that local gates are technologically easy compared to non-local 
gates, the curious reader may wonder why bother with distributed quantum com- 
puting? Surely we could build a quantum computer using these techniques for local 
gates? Unfortunately this is often not the case. Arbitrary control over a few qubits 
at each node is much more technologically viable than arbitrary control of an ar- 
bitrary number of qubits. If we try and scale up nodes to larger sizes, numerous 
technological difficulties emerge. One of these is the problem of addressability in 
completely local quantum computers, w hich occurs when we only want one qubit 
to interact with a particular control field-*-* 1 . This can not be performed by spatial 
resolution as interacting qubit systems will typically be separated by less than an 
optical wavelength. As for frequency resolution, this may be possible for a limited 
number of qubits. For example, in a 2 qubit system frequency splitting may oc- 
cur because the qubits are of different species, such as electron and nuclear spin. 
However, in long chains of identical qubits there will be no effect that splits the 
frequency of every qubit by a resolvable amount. 

Although addressability is an obstacle to large-scale local architectures, it is 
certainly not an insurmountable obstacle. For example, there exist so-called global 
control schemes that circumvent the addressability 

pr blem^l^|76|77|78| 

Even 

when addressability can be achieved, cooling presents an additional hurdle, as it 
is often extremely difficult to reduce population of low lying excited states in such 
systems. 

Aside from the question of technological viability, distributed quantum com- 
puters have other inherent advantages. Typically, a distributed archi tectur e can 
be made fault tolerant at higher noise levels than a local architect ure eP^ QI p ur _ 
thermore, distributed architectures allowing for arbitr ary connectivity offer com- 
putational performance superior to local architectures^. As such, it may prove 
advantageous to develop a distributed architecture even if a local architecture is 
viable. Having outlined some issues involved in local operations, and local architec- 
tures, we return to our main concern: how to generate entanglement between local 
nodes. 

3.4. Single-photon protocols 

Cabrillo et aff^l proposed the first protocol for generating entanglement between 
distributed atoms via photon measurements, and this seminal work has spawned nu- 
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Fig. 8. single-photon protocols: (a) The experimental set-up used for single-photon protocols 
(and many other protocols). Two matter systems embodying separated qubits are placed so that 
emitted photons are encouraged to emit towards a beam splitter, this can be achieved by using 
QED cavities or collection lens. Before the beam splitter are frequency filters that absorb light 
emitted from the |e) — * |0) transition. After the beam splitter are two photon detectors, (b) The 
A-level structure required for each matter system in the original Cabrillo proposal. 

merous spin-off protocols. Bose et a P showed that this protocol could be adapted 
for the purpose of implementing teleportation. Browne et a P used adiabatic laser 
pulses to improve the protocol's success probability and robustness to various er- 
rortQ. Childress et demonstrated that a different electronic level structure can 
provide the same benefits. Here we will outline the essential components of the 
Cabrillo proposal. For the purpose of better comparison with later proposals we 
will assume that two photon detectors are use 

Cabrillo et al considered two matter systems with a A-level structure, with an 
excited state |e) that decays to ground states |0) and |1) with decay rates T and Ti 
respectively (see Fig. [5Jd). These matter systems are arranged in the experimental 
set-up shown in Fig. [8^,. Note that frequency filters ensure that only photons from 
the |e) — > |1) can reach the detectors. The Cabrillo proposal begins with both 
matter systems prepared in the state |0). Next a short pulse of laser light resonant 
with the |0) «-> |e) transition is used to pump a small proportion of the population 
into the state |e), giving the joint state: 



where 9 is a small parameter that determines the strength of the weak excitation. 
Since there is only a small probability of excitation, usually no photons will be de- 
tected. However, when a photon is detected we project out the 10)^ |0) B component. 
The purpose of choosing small 8 is so we can assume that two photon processes 
are negligible, so \e) A \e) B also vanishes. This is necessary because although only 
one photon was detected, photon loss means that more could have been emitted. 

'For example, by using detuned lasers they can filter out scattered laser light. 

m In the original proposal only a single photon detector was proposed. This wastes an output port 

of the beam splitter and reduces the success probability by 1/2. 



(cos(0) |0) A + sin(0) |e) A )(cos(0) |0) fl + sin(0) \e) B ), 



(40) 
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Later we derive the infidelity due to non-zero though for now we assume the two 
photon term is negligible. The presence of the beam splitter ensures that we do not 
know where the emission originated, and hence we are in a maximally entangled 
Bell state: 

|**) = -L(|0)|l) + e < *|l)|0)) > (41) 

where the variable <fi depends on the path lengths to the beam splitter and which 
detector clicked. All single photon schemes require that the apparatus is interfcro- 
metrically stable to distances much shorter than an optical wavelength so that the 
phase factor <f> is known. 

Now we account for the infidelity due to non-zero 9, which corresponds to double 
excitation processes. If two photons are emitted, and one is undetected, it will 
eventually be absorbed by the environment. Equivalcntly, both may reach a given 
detector which is unable to resolve, or count, the pair. Either process is described 
by a measurement of the number of photons, but without our knowledge of whether 
the answer was 1 or 2. Hence, we have an incoherent mixture of the desired state 
with 1 11): 

(1 - V )\*+)(*+\ + 77111X111, (42) 

where rj depends on the excitation strength and photon loss rate, and for simplicity 
we have assumed that T <C Ti. The ratio of the weighting (1 — 77) : 77 is proportional 
to the relative probability of one and two photons resulting in a single detector 
click. Defining T as the probability of transmission and detection of a given photon 
emitted by either matter system, it follows that one photon emission and detection 
occurs with probability: 

P (1) = 2 sin 2 (6>) cos 2 (9)T. (43) 

The equivalent probability for two photon emission followed by one detector click 
is: 

P (a) =sin 4 (0)(l-(1-T) a ), (44) 

where we have assumed that the photons bunch and detectors can not resolve the 
number of photons. It follows that the erroneous contribution has magnitude: 

P (2) =ri n»(0)(2-D 
P(i)+P(2) 2- sin 2 (9)T V ' 

To reduce rj we must reduce 8, and hence the probability of detecting a photon. In 
the small 9 and small T regime, we have direct proportionally between infidelity 
and success probability: 77 ~ 9 2 and Pm + P( 2 ) ~ 26* 2 . 

Note that in Cabrillo et aVs original proposal the decay channel To only reduces 
the success probability. Therefore, we can also use a similar approach to entangle- 
ment generation where only the state |1) couples to |e). For such a scheme we 
prepare the superposition cos(6>) |0) + sin(6>) |1) and then ir pump |1) — * |e). 
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Entanglement from a single p hoton detection has been demonstrated between 
two clouds of cold caesium atomdM122l j n ^ese experiments the logical states of 
the qubits were the: |0) all atoms in ground state; |1) an ensemble of atoms sharing 
an excitation. 



3.5. Two photon protocols 

Unlike single-photon schemes, the later two-photon schemes are robust against 
interferometric instability and do not have their success probability tied to infi- 
delity. As outlined earlier we will introduce a typical anti-bunching protocol and 
the double-heralding protocol. It is a testament to the promise of distributed quan- 
tum computing that an anti-bunching scheme h as succ essfully been used to generate 
entanglement between two remote trapped i on iP^j ^ . 

The repeat-until-success (RUS) protoco ff^^ is another interesting protocol 
that we do not have space to review. However, it is worth noting that the RUS pro- 
tocol performs a generalized measurement with outcomes that arc locally equiv- 
alent to the identity or control-^ gates. Therefore, it is the only known hybrid 
protocol that has no failure^] outcome when there is no photon loss and detectors 
can resolve different numbers of photons. 



3.5.1. Bunching based schemes 

Two identical photons impinging on different different ports of a beam splitter will 
bunch together and leave through the same beam splitter port (see section 13. ip . 
This feature of Bell states has been ex ploited in numerous proposals for generating 
entanglement either using polarizatio n 1 ^^ or frequencjE^ to define the logical 
basis. All of these schemes essentially use a A-level structure and success is heralded 
by anti-bunched photon detections, from which we can infer that the phot ons were 
not identical. Here we will describe the proposal of Duan and Kimble^!. 

Duan and Kimble proposed an experimental set-up that uses waveplates and 
polarizing beam splitters to erase which-path information of polarization encoded 
photons (see Fig. [9]). Their matter system has a level structure with an excited 
state that decays to |0) or |1) via the emission for a horizontal or vertical photon, 
respectively. They also assume a ground state \g) that can be driven to the state 
\e). The role of \g) is only as a storage state before driving the state to |e), so it 
is unnecessary if we have some other means of preparing |e). The protocol runs as 
follows: 

n By generalized measurement we mean Positive Operator Values Measurement, or POVM for 
short 

°The identity operation is not a successful entangling gate. However, if the qubits hold some 
pre-existing entanglement with third parties, this third-party entanglement is not destroyed by 
the identity gate. In this context, the identity gate can be interpreted as neither a success nor a 
failure. 

p We give their more refined version of the protocol that uses more experimental apparatus but 
increases the probability of success. 
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waveplate(7r/2, 0) 

Fig. 9. An outline of Duan and Kimble's proposal for generating entanglement that is heralded by 
detecting anti-bunched photons, (a) the experimental apparatus includes 3 waveplates, 3 polarizing 
beam splitters (PBS) and four photon detectors. Dli an( i Dl2 are tire "left" detectors and the 
"right" detectors are Dri and Dr2- (b) the electronic level structure. 

(1) From the initial state \g) \g), adiabatically drive both matter-cavity systems 
into the state |e) \e); 

(2) wait for the system to relax whilst monitoring for emitted photons; the protocol 
succeeds when we observe two photon detection events, one in a left detector 
and one in a right detector. 

Just after the matter systems A and B relax and emit a photon, the state of the 
system is: 

H') = (go |0> at + 9l |1) at)(g |0) b\ + 9l |1) b\) |vac) , (46) 

h v h v 

where go and 171 are the relative probabilities of relaxation by horizontal versus 
vertically polarized photons. The first waveplate rotated the polarization at mode 
b: 

|V) = (90 |0) 4 +5i |1) ot)(ff |0) b\ - 9l |1) bU |vac) . (47) 

h u u h 

Impacting on the PBS the vertical photons are reflected: 

H'} = ("3o |00> clcl - .90.91 |01) cUl + 3130 |10> 44 + A I 11 ) l vac > - ( 48 ) 

h v h h v v h 

for which the subspace with one photon on either side of the PBS is entangled: 

l^>anti-bu„cncd « 3i3o(|10) M - |01) cU) |vac) . (49) 

However, when we measure the photons we have to ensure that neither we nor 
any other system has information about whether the photons are h or v polarized. 
If this information is measured by any system then it will result in a projection 
into 1 01) or |10). For this reason we use another waveplate and PBS to measure 
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modes c and d in the 45° polarization basis {(^r), (^/^)}> which gives completely 

random outcomes for h or v photons. Hence, a successful application of the Duan 
and Kimble protocol will project the matter system into a maximally entangled 
statci. 

Compared to single photon protocols, antibunching schemes are more sensitive 
to photon loss as they require two photons to successfully reach the detectors and 
may have a more complex apparatus and hence a higher loss rate. However, photon 
loss only affects the success probability and not the fidelity. Furthermore, since 
each term of the superposition contains a photon that traverses both arms of the 
apparatus, the scheme is insensitive to differences in path lengths. 

3.5.2. The double-heralding scheme 

We now turn to the double-heralding scheme^! which uses matter qubits with an 
L-level structure arranged as shown in Fig. 1101 The procedure consists of four steps: 
preparation, first heralding, bit flips and second heralding: 

(1) Preparation: Both qubits are initialized in the |+) = (|0) + 11))/-^ state; 

(2) First heralding: Both qubits are subjected to a 7r-pulsc of frequency lu, which 
causes the transition |1) — > |e). Next, the photon detectors D\ and Do wait 
for detector clicks until enough time has passed that it is very unlikely that 
the cavity still contains any excitations. We proceed to the next step of the 
protocol only if we have seen exactly one click. 

(3) Bit flips: Both qubits are bit flipped by the rotation XaXb; 

(4) Second heralding: We repeat the 2 nd step. If and only if we again see one click 
then the protocol has succeeded. 

Notice that the first two steps are the same as for the Cabrillo protocol when 
9 = 7r/4 and decay is only possible by one channel. Hence, after the first heralding 
we have a mixed state: 

(db) l * 0><VI/ ' 1 + (f=?) |11>(n| - (50) 

After the bit flip we have 

(dbO i^k^i + (fr-f ) io°><o°i- (51) 

When we perform the second heralding, the erroneous term can not emit any pho- 
tons and so can not lead to a detector click. Furthermore, on the second heralding 
the photon traverses the opposite path to that of its first journey, and hence both 
terms of the superposition acquire the same phase due to path length difference^]. 

q There may be a trivial phase depending on which detectors click. 

r This assumes that path length has not varied between the first and second round. 
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Fig. 10. An outline of Barrett and Kok's protocol for generating entanglement by double herald- 
ing the parity measurement made by a single photon scheme, (a) the experimental apparatus is 
similar to that in the Cabrillo proposal, shown in Fig. [8] but with the removal of the frequency 
filters, (b) an L-level electronic structure; |0) and |1) are degenerate, or nearly degenerate, in en- 
ergy; The excited state |e) decays to |1) via a photon emission, whereas decay /absorption between 
|0) <-> |e) is not permitted due to a selection rule. 

Although, in any one successful round only a single photon is detected, the whole 
protocol is successful after two photons have been detected. 

Double heralding is not just suitable for producing Bell pairs, but can also be 
used to produce a parity projection for growing graph states (see section FfTTj) . For 
graph state growth we omit the preparation step, and apply the protocol directly 
to the graph state qubits. Upon success we undo the XaXb, and the result is an 
odd parity projection. 

3.6. A scaled up device 

To recap, in this section we have reviewed some protocols for generating entangle- 
ment and the physical mechanisms they use. Although entanglement is essential 
for quantum computing, quantum states with bipartite entanglement alone are not 
sufficient. Rather we require that many nodes of a distributed quantum computer 
can be entangled with each other. Therefore, all of the entanglement proposals will 
require an optical multiplexer that switches the paths of photons, as illustrated in 
Fig. nu 

In the case of double heralding we are able to grow graph states by omitting 
the preparation step. When successful, the double heralding procedure produces an 
entangling parity projection, but when it fails it projects two qubits into a separable 
states. These success and failure outcomes have a simple graphical interpretation in 
terms of either connecting graph vertices or removing two vertices from the graph. 
Even though double heralding grows graphs probabilistically, we see in section |4~T1 
that strategies exist for efficient graph growth. 

However, what holds for double heralding will not automatically transfer over to 
other schemes for generating entanglement. For example, A level schemes require us 
to initialise in some separable state that can decay to either logical state, and this 
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Fig. 11. Schematic of a DQC: Nodes of a DQC are alphabetically labelled, and each node contains 
some number of qubits of which at least one is optically active. Within a node control fields can 
manipulate qubits to high fidelity, so each node is miniature quantum computer. However, large 
scale quantum computing requires entanglement generation between nodes. Photon decay paths 
from different nodes are switched by an optical router, enabling entanglement generation between 
any two nodes. 

loss of prior entanglement during the preparation is unavoidable. Thus, to make 
use of such protocols it is necessary to have at least one additional qubit present 
at each node, where the two (or more) qubits within a node interact in such a way 
that we can perform good quality quantum gates between them. Then the optically 
active qubit s can neg otiate entanglement, while the secondary qubits act as the 
logical It is well known that once a Bell pair is shared between two 

nodes it can be used to implement an entangling two-qubit gate, and hence enable 
either circuit model quantum computing or one-way quantum computing. When 
using ancillary qubits in this fashion the graph state is also immune from damage, 
reducing the need for strategies for probabilistic graph growth, but at the cost of 
additional node structure. 



4. Growth strategies 

4.1. Graph growth with parity projections 

From the constructive definition of a graph state, it trivially follows that if we 
can perform CZ gates and local operations then we can construct graph states. 
However, many proposals for entanglement generation are projective rather than 
unitary, especially when the entanglement is ind uce d by measurements. Indeed, to 
our knowledge the Repeat-Until-Success protocoff^is the only measurement based 
protocol for generating entanglement that upon success generates a true control- 
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Z gate as opposed to a projection. Whereas most of the entanglement protocols 
for distributed quantum computing only succeed when we detect some odd parity 
photon measurement signature, which we might expect corresponds to an odd parity 
projection on some qubits x and y: 

P- = |01)(01| 1BiS + |10)(10| x , 1 , i (52) 

where we use the negative sign to denote odd parity because it can also be written 
in term of Pauli-opcrators: 

P ± = {1±Z X Z V )/2 1 (53) 

which also defines the even parity projector. P + . Since ±Z x Z y is a stabilizing 
operator, this notation shows that P ± project onto the subspace stabilized by the 
±Z x Z y of the same sign. Since such parity projectors project onto a stabilizer state 
we might expect that they are suitable for constructing graph states. Indeed, in 
Fig. [12k we illustrate how a parity projection affects two qubits A and B: they 
become redundantly encoded, which is LU-equivalent to a graph state with qubit A 
gaining all of the neighbours of qubit B, whilst B becomes a dangling cherry with 
A as its only neighbour. 

To verify that this is the effect of an odd parity projection, we first consider the 
initial state of the constructive defined graph state, which we divide into three 
partitions: qubit A, qubit B, and the remaining graph state |Sf'): 

\%) = |(|0) + |1) J^M|0) + |1) 2f B ) B |Sf') , (54) 

where 3f x denotes a product of all the phase-flip operators resulting from the ap- 
plication of control- Z gates involving graph vertex x, or formally: 

It follows that projection into the odd parity subspace generates: 

|Sf/) oc P~ |Sf<) , (56) 
= i=(|0,l}#B + |l,0> 3T a )a,b |Sf'). 

= -j=($? B X B H B )(\0, +) + \l, -) % A 2?b)a,b |Sf'> , 

= ±=(!¥ B X B H B ) (CZ£(|0) + |1) 2? a 2?b)a \+) b ) \&') , 

where removal of the local unitaries 3f B X B H B clearly gives the constructively de- 
fined graph state depicted in Fig. WZk . Note that we have only indicated the local 
unitary H B in the figure. Throughout this review we will omit Pauli-opcrators from 
our illustration of states LU equivalent to some graph state. We use this convention 
because the presence of a Pauli-opcrator does not dramatically effect the evolution 
of a graph state under Pauli measurements, whereas the presence of a Hadamard 
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Fig. 12. Examples of graph growth by parity projection: (a) the outcome of a parity projection 
on two qubits both part of a graph state with many neighbouring vertices. One qubit, say A, 
acquires all of the neighbours of the other qubit, B. Whilst qubit B becomes connected to qubit 
A and only qubit A, with the addition of a Hadamard rotation. We can remove this Hadamard to 
obtain a constructively defined graph state. When a qubit only has 1 neighbour, as qubit B does, 
we call this qubit a cherry. Note that parity projections do not have a basis towards one qubit, 
we could equally describe the resulting state with qubit A as the cherry. Since the roles of A and 
B are interchangeable they are said to be redundantly encoded, (b) the growth a chain of qubits 
with a probabilistic control- Z operation. Provided an offline resource of chains greater than some 
critical length, the longest available chain will, on average, grow in length. 

rotation does have a significant effeclQ. Note that an even parity projection would 
give the same state up to a difference in local Pauli operators. 



4.2. Probabilistic graph growth 

A broad range of entanglement protocols, generating control- Z gates or parity pro- 
jections, enable graph growth. However, entanglement generation in distributed 
quantum computers will typically be probabilistic, with a failed attempt projecting 
the qubits onto some separable state. This raises the question of whether proba- 
bilistic graph state growth can be performed efficiently. A naive approach to graph 
growth would restart whenever it fails to entangle any two qubits. Denoting the 
probability of generating entanglement by p, the overall probability of success for 
this naive strategy will fall exponentially as p n , where n is the number of required 
number of attempts at entanglement generation. An exponential fall in success 

s For example, a graph state qubit with a Hadamard measured in the X-basis will generate the same 
final state as a Z-basis measurement on an equivalent graph state without the Hadamard present. 
In contrast Pauli-operators are only capable of more interchanging the outcomes of Clifford group 
measurements . 
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probability entails an exponential amount of time until success, and hence we lose 
any speed-up gained from using a quantum algorithm. Therefore, we say that a 
strategy for probabilistic graph growth is efficient when the temporal and spatial 
overheads scale polynomially in the size of the graph state. 

Indeed numerous efficient strategics do exist, and they all explo it the way failed 
entanglement generation causes only local damage to graph state d I I I I Typi- 
cally, the exposition of these strategies has assumed that entanglement is produced 
by a control- Z gate, though it is trivial to modify them for graph growth by prob- 
abilistic parity projection. For the entanglement generation protocols that we are 
interested in, failure to entangle two qubits will always either: project the qubits into 
a separable state of the computational basis or produce an unknown phase error 
on the qubits. Therefore, we can measure out these qubits in the Z-basis, applying 
a correction to their neighbours where necessarj@. We say that failure damages the 
two qubits, and that they have to be reset before we attempt to rcgrow the lost 
graph vertices. Note that if a failure could cause bit-flip noise, we would also have 
to reset their neighbourly. This property of graph states is closely related to a trick 
introduced into linear optical quantum computing by Knill et aP ^j, w hich teleports 
successful instances of probabilistic gates into the main algorithm^. 



4.2.1. Dynamic growth strategies 

Here we outline two strategies that we call dynamic graph state growth strategies, 
so-called because they do not use a predecided sequence of attempts at generating 
entanglement. Implementing these protocols would require our quantum computer 
to be able to attempt entangling operations between any, or almost any, pair of 
qubits. First we consider the growth of long chains of qubits, which is the first 
stage of graph growth for the proposal of Duan and RaussendorP. In Fig. If 2b we 
show the failure and outcome possibilities when connecting two chains, each initially 
containing L qubits. This succeeds with probability p, and the successfully joined 
chain has length 2L. However, if it fails, then the result is two chains of length 
L — 1. From this, we can calculate the expectation value of the longest available 
chain: {Li) = p.2L + (I — p)(L — I). For, on average, growth of the chain, the 
expected final chain length must be larger than the initial chain, (L) < (Lf). This 
simplifies to requiring that the initial chains are longer than some critical length 
L c = p^ 1 — I. Growing these critical length chains can be expensive for small p, but 
the resource cost does not scale with the size of the computation. Therefore chain 
preparation is a constant cost that does not affect the efficiently of the strategy, and 

'A |0) result requires no correction, whereas a |1) outcome requires that that all graph neighbours 
are phase-flipped. 

"More formally, if qubit x is involved in a failed procedure resulting in bit-flip noise then we would 
have to reset all qubits y such that y S N(x). Duan and Raussendorf assumed this worse case 
scenario when evaluating their graph growth strategy^]. Therefore, when we describe their strategy 
in the next section there will be slight differences in some expressions 
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Fig. 13. Examples of Re-routing growth strategies: (a) the cross strategy proposed by Duan and 
Raussendorf: (b) the micro- cluster strategy proposed by Nielsen. For an idea of the number of 
qubits required, the number of qubits/connections lost at each stage of graph growth has for a 
success probability of p = 1/3. 



we call such states offline resources. Furthermore, Duan and Raussendorf showed 
the time to prepare chains of length L scales as log 2 (i — L c ). From here Duan and 
Raussendorf propose that chains arc connected at the middle to produce a large 
cross shaped graph state, as outlined in Fig. [T3k . The central vertex of this cross 
will be present in the final target graph state, the purpose of all other qubits is 
simply to buffer against graph damage. When attempting to connect two crosses 
the probability of an outright failure drops exponentially with the size of the buffer. 
After a successful connection, any remaining buffer qubits are removed by Clifford 
group measurements in the X or Y basis. By constructing suitable sized cross 
shaped graphs states, the chance of failing any connection in building an N by N 
square lattice, or cluster state, is kept to a small probability at subcxponential cost 
in time and qubits^. A strategy in a similar spirit is the micro-cluster strategy 
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Fig. 14. A simple example of graph state growth by percolation. First we produce microclusters 
arranged in a grid, and then attempt to join the microclusters into a square grid, much like the 
strategies in Fig. 113b . However, rather than ensuring that the microclusters are large enough 
that every connection is almost certain to succeed, many will fail. With the grid divided into 
blocks of qubits (we show 4 blocks with different colours which slightly overlap), we require that 
a qubit can be found in each block such that a path of successful connects leads from that qubit 
to neighbouring qubits. When the successful probability is greater than some threshold (1/2 for 
square lattice and 5 qubit microclusters), it is asymptotically certain that this strategy succeeds. 



proposed by NielserP, which wc outline in Fig. [T3b . Instead of cross shaped graph 
states, this strategy proposes star shaped objects known as micro-clusters, of which 
only the central vertex will be present in the final graph state. This strategy has a 
similar scaling to that of Duan and Raussendorf. Note that micro-cluster are LU- 
equivalent to GHZ states, named after Grecnbcrgcr, Home and Zollinger who first 
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described them^ 1 * and as such the terminology often used interchangeably. 



4.2.2. Percolation strategies 

Both the aforementioned strategies assume that the quantum computer is capa- 
ble of attempting an entanglement operation between almost any two qubits, with 
the decision making process made dynamically by an assisting classical computer. 
Building a quantum computer that provides the flexibility of being able to change 
which qubits we entangle is clearly more technologically challenging than one that 
does not. A natural question is whether graph states can be grown when the quan- 
tum computer can only entangle each qubit with a fixed number of prcdecided 
qubits? Indeed it has been shown that percolation phenomena can b e exploi ted 
to grow graph states when the intcrconncctivity of a device is limited^ lO I ^M In 
Fig. [14] wc illustrate a simple percolation strategy where no qubit needs to be en- 
tangled with more than 4 other qubits. The overall square lattice is sub-divided 
into partially overlapping square blocks that will each represent one qubit in the 
final graph state. Each block is successful if there is a uninterrupted path from left- 



to-right and top-to-bottom, and we also require crossings in the overlapping region. 
Assuming that such paths exist, all qubits except one in each block are measured in 
a Pauli-basis to reduce down to a square lattice. The percolation phenomena that 
they exploit is well known in other fields and describes phase changes in problems 
such as the movement of solvents through filter paper (chromatography) and the 
movement of petroleum through fractured rock. Kicling et ai^ utilise results from 
bond-percolation theory that to show that provided the probability of a hole de- 
fect is smaller than some critical value, then at a polynomial overhead the success 
probability is asymptotically certairQ. For a square lattice this critical probability 
is 1/2, though more favourable probabilities can be tolerated by utilising different 
lattice topologies, such as 3D square lattices and diamond lattices, and by simply 
increasing the size of the initial microclusters. Perc olation phenomena has also been 
investigated when sites rather than bonds faiff^, and in the context of quantum 
communication^! 



4.2.3. Problems with probabilistic graph growth 

All the above strategies take the probability of successful entanglement generation, 
p, as a constant, and assume perfect fidelity operations. Given these constraints it 
is claimed that graph growth can be performed efficiently. However, as p decreases 
the overhead in time and qubits required increases faster than I /p. Indeed, ex- 
tremely low p makes graph growth impossible when the expected time to make a 
single successful connection exceeds the decoherence timfQ of the qubits. Paying 
heed to the effects of noise, as we increase the number of qubits used in a single 
cross, micro-cluster or percolation block, we i ncre ase the amount of noise that will 
eventually collect on a single graph state qubitl^H Hence, under realistic conditions 
probabilistic graph growth strategics will be too noisy for quantum computing when 
p is too low. Exact bounds are not known, but what is clear is that significant gains 
can be reaped by generating entanglement with a lower probability of damaging 
the underlying graph. 



5. Robust computation 

While the growth strategies discussed in the previous section allow us to deal with 
probabilistic entangling operations, more effort is necessary to deal with noise ac- 
cumulated during construction and measurement of graph states. Such errors may 
be caused by errors in the local operations or entangling operation, for example due 
to photon detector dark counts, or simply by decoherence of the matter qubits used 
to store the graph state. In any case, any scalable architecture must take steps to 



v If the target graph state is a square lattice of L by L qubits, and we use blocks of order L M by 

L M qubits ( fj, > 0), then the probability of success approaches 1 as L — > oo. 

w This is a measure of the rate at which noise takes a qubit into a randomized state. 
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suppress such noise in order to maintain a finite probability of obtaining the correct 
output from a computation. 

5.1. Error correction and fault- tolerance 

One of the major successes of quantum information theory has been the devel- 
opment of a formal theory of fault-tolerant computation. In such a model, the 
quantum information is encoded using an error correction code to protect against 
interaction with the environment and other random errors. A universal set of gates 
is then constructed in such a way that errors are not amplified. A regular regime of 
error correction between logic gates ensures that the computation remains in the 
encoded subspace. Provided that the error rate is below a critical value, known as 
the fault-tolerance threshold, errors can be suppressed indefinitely with only poly- 
nomial resource overhead. The exact value of the fault-tolerance threshold depends 
on several factors including the specific error correction code used, the structure of 
the fault-tolerant logic gates and on the underlying physical operations available. 

Maximising the fault-tolerance threshold is an extremely important open prob- 
lem as it allows for computation to be performed in noisier systems, reducing 
the technological overhead required to implement quantum computation. Although 
many families of quantum error correcting codes have been proposed, determining 
a threshold is far from trivial. A number of different factors affect the threshold 
for fault-tolerance. Important factors include how many errors can be detected and 
diagnosed by an error-correcting code and how many gates are required per round 
of error-correction. ^^^^ 

While standard error correction codes, such as the Shor and Steane codes™, 
can be applied to one-way quantum computing, it is far from clear whether codes 
designed for the circuit model of quant um computation will prove to be the best 
option in this new ar 

One promising proposal, brought forward by Raussendori*^, is to use surface 
codeJ^ to protect against errors. Surface codes represent a topological approach 
to error correction, in which a tiling of local operators on some surface are used 
to stabilize the subspace in which the computation takes place. An advantage of 
this approach is that it finds a natural substrate in three dimensional extension of 
cluster states. Such topological approaches to fault tolerance provide a fertile area 
for further research and may well result in higher fault tolerance thresholds. 

5.2. Entanglement distillation 

We have briefly reviewed how fault tolerant quantum computation enables reliable 
computations despite the threat of errors occurring at any stage, e.g. graph state 
growth, storage, or measurement. However, these techniques can only cope with 
small error rates. But what if noise is significantly worse at a particular stage, such 
as during our entangling operation used for graph state growth. This possibility 
naturally occurs in a distributed architecture, where entangling operations between 
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nodes will occur via a fundamentally different mechanism than inside nodes. For hy- 
brid matter/optical systems the entangling operation may be noisy dues to photon 
loss, dark counts^), imperfect alignment of optical elements or non-identical photon 
sources. Whilst there will also be noise processes within nodes, it seems that these 
must strictly be less noisy than the appropriate fault tolerance threshold. Provided 
that noise is low inside nodes, many uses of a noisy inter-node entangling operation 
can be used to produce a single high fidelity entangling operation. The family of 
protocols that achieve this arc aptly called entanglement distillation protocols, and 
can provide high fidelity graph state growth even when the entangling operation is 
significantly nosier than any fault-tolerance threshold. 

Such a powerful tool does not come for free and requires some additional com- 
plexity in the structure of quantum computer. To employ entanglement distillation 
protocols requires more than one qubit per node, with the precise number of qubits 
depending on the type of noise and protocol employed. Most of the literature focuses 
on so-called Pauli noise, where the entangling operations prod uces a Bell state up 
to random Pauli Z (phase), X (bit) or Y errors. Dur and BriegeP^ showed that 4-5 
qubits per node can be used to dist ill a lmost any entangling operation that suffers 
from Pauli noise. Later Jiang et a|E2 proposed a different scheme that distilled 
general Pauli noise using 4 qubits per node, but only needed 3 qubits per node 
for phase noise. Both these protocols can be employed both within the one-way 
model and the circuit model, as they use a one qubit buffer per node to store a 
high fidelity Bell pair which is then consumed to perform any high fidelity logical 
operation. Exploiting the simplicity of graph state growth, it has been showr r" 1 1 
that the need for a one qubit buffer can be eliminated, with high fidelity graph 
states growth possible using only 2 (phase noise or photon loss noise) or 3 (gen- 
eral Pauli noise) qubits per node. A description of the internal working of these 
protocols is beyond our current scope, but we hope to have conveyed that a great 
deal of inter-node noise can be distilled with a node containing a modest number of 
qubits. Note that, if much larger no des are a vailable then we may use more complex 
protocols for graph state distillatiorl22UflQI. 



6. Conclusions 

In this introduction wc have seen how hybrid matter-optical systems can be com- 
bined with the one-way model of computation to yield an extremely powerful and 
robust architecture for scalable quantum computing. As such technologies mature, 
additional local structure will allow for higher efficiency, by buffering growth and 
distilling the entanglement present. These systems hold the promise of universal 
control without the limitations imposed by local architectures, and may offer a 
shortcut on the route to scalable devices. 



x Where photon detectors click in absence of a photon. 
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